Filling the void in confined polymer nematics: phase transitions in a minimal model 

of dsDNA packing 



O 



(N 
(N 



Homin Shin and Gregory M. Grason 
Department of Polymer Science and Engineering, 
University of Massachusetts, Amherst, MA 01003 

Inspired to understand the complex spectrum of space-filling organizations the dsDNA genome 
within the capsid of bacterial viruses, we study a minimal, coarse-grained model of single chains 
densely-packed into a finite spherical volume. We build the three basic elements of the model- 
i) the absence of chain ends ii) the tendency of parallel-strand alignment and iii) a preference of 
uniform areal density of chain segments-into a polymer nematic theory for confined chains. Given 
the geometric constraints of the problem, we show that axially symmetric packings fall into one of 
three topologies: the coaxial spool; the simple solenoid; and the twisted-solenoid. Among these, 
only the twisted-solenoid fills the volume without the presence of line- like disclinations, or voids, and 
are therefore generically preferred in the incompressible limit. An analysis of the thermodynamics 
behavior of this simple model reveals a rich behavior, a generic sequence of phases from the empty 
state for sinall container sizes, to the coaxial spool configuration at intermediate sizes, ultimately 
giving way, via a second-order, symmetry-breaking transition, to the twisted-solenoid structure 
above a critical sphere size. 
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I. INTRODUCTION 



Consisting primarily of a protective capsid enclosing 
a genome that codes for minimal number of proteins, 
viruses are among the simplest living organisms. Given 
the economy of components to the virus, it is remarkable 
that successful replication relies on the robust solution of 
a rather complex materials geometry problem: packing of 
a long chain molecule within an enclosed volume of linear 
dimension that is nearly 1000-fold smaller than the chain 
length. This problem is particularly pronounced for bac- 
terial viruses which carry a dsDNA genome, whose per- 
sistence length, 50 nm, is comparable in size to the inner 
radius of many bacteriophage capsids [l| . Understanding 
the structure, kinetics and mechanics of is not only of 
importance for understanding this crucial stage the viral 
life cycle, it also poses fundamental challenges regarding 
the underlying frustration of packing a one-dimensional 
object within a finite, three-dimensional volume at max- 
imum density. 

To account for the huge pressures (~ 10 atm) needed to 
confine a dsDNA genome within the viral capsid, a num- 
ber of theoretical studies [2|-|3| have invoked an idealized 
"spool" configuration [^ in which the chain winds axially 
around a central and largely empty central core. While 
this configuration is convenient for quantitative study of 
the energetics of DNA packing, a number of recent the- 
oretical and experimental studies suggest the spacing- 
filling structure of DNA within capsids is often far more 
complex. Simulations of coarse-grained chains forced into 
empty icosohedral or spherical volumes [10l4l3j find com- 
plex chain packings with markedly lower symmetry than 
the spool, twisted or folded toroidal structures [l4|. Most 
striking are the results of high-resolution EM studies 
of the fully encapsidated T5 bacteriophage [l5[, which 
find that DNA is present at high density and locally 
well ordered throughout the volume of the capsid, but 



modulated in a complex global pattern to accommodate 
the constraints of confinement. Despite intense inter- 
est devoted to DNA packing in recent years, a number 
of elementary questions about the structure of confined, 
densely-packed chains remain unanswered. What, pre- 
cisely, is the underlying source of frustration, what struc- 
tures best negotiate the competing geometrical interests 
of the problem, and what are the microscopic and ther- 
modynamical parameters that discriminate among these? 
In this Letter, we describe an approach to this problem 
based on a minimal description of the frustrated pack- 
ing of a single chain molecule into a spherical volume: 
the coarse-grained theory of polymer nematics. We pro- 
pose a model consisting of the following three elements 
i) a vanishing density of chain ends; ii) a preference for 
parallel orientation of neighbor strands; and iii) a pre- 
ferred non-zero areal density of chains throughout the 
volume. The first of these ingredients, a geometrical con- 
straint, is described in terms of conserved, or divergence 
free, order parameter, the chain fiux t(r), whose direction 
and magnitude, respectively, encode the mean values of 
the chain orientation and the local density of chains in 
a perpendicular plane. Based on this constraint alone, 
we classify all axially symmetric packings into one of 
three distinct topologies. Among these, only a "twisted- 
solenoidal" structure achieves a non-zero flux throughout 
the volume, while the coaxial spool and simple solenoid 
structures require the presence of line-like voids, or discli- 
nations. Using variational and numerical methods, we 
study the simplest possible phenomcnological free en- 
ergy that describes the thermodynamic competition be- 
tween different states of packing. While coarse-grained, 
polymer nematic models have been applied to this prob- 
lem before 0, [l6l - [l8j . the key and novel result of this 
study is the complete determination of the phase behav- 
ior in terms of a single parameter, the container radius R, 
which quantifies the relative preference for uniform pack- 
ing over the elastic cost of non-uniform chain orientation. 




chain flux must satisfy tfic following conservation law, 
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FIG. 1. The 3 topologies of spherically-confined, azimuthally 
symmetric, end-free chain configurations: (a) coaxial spool; 
(b) simple solenoid; and (c) twisted solenoid. Streamlines 
denoting the orientation of chain flux are visualized on nested 
toroidal surfaces of constant, ^(r), the Stokes stream function 
of the solenoidal configurations. While lines of t(r) in (a) 
and (b) become singular along line-like disclinations-along the 
axis of the spool and along the circular center of the solenoid, 
respectively-the flux conflguration of (c) is singular at only 
at point-like defects on the poles of the spherical boundary. 



We find a universal sequence of equilibrium phases with 
increasing R: the empty — >■ spool — > twisted solenoid. Fi- 
nally, we explore the a higher-order theory that accounts 
for the distinct costs of twist and bend deformations of 
the chain orientation, leading to stability of the simple 
solenoid packing for intermediate R. 



II. GEOMETRY OF END-FREE CHAIN 
CONFIGURATIONS 

To describe configurations of single chains confined 
within a spherical volume, we introduce a coarse-grained 
vector field corresponding to the local polymer "current" , 
t(r) — p(r)n(r). Here, n(r) is the mean orientation or 
chains segments at r and p{r) is the mean flux of chains 
in a plane perpendicular to n(r). It is important to 
note that p{r) carries information both about mean areal 
density of chains, as well as the degree of orientational 
order at r. In principle, this field may be further de- 
composed as, p{r) = pm(r)s(r), where Pm{i^) is the local 
area-densit y o f chains and s(r) is the local vector order 
parameter |18| . 

In densely-packed configurations of a single chain, ends 
are rare. Here, we assume that the idealized geometry 
of dense, single-chain packings is characterized by the 
extreme limit where the density of chain ends vanishes 
throughout the volume. In the absence of chain ends, the 



V-t = 0. 



(1) 



At the boundary of the volume, at r = i?, the absence of 
chain ends also requires the vanishing of outward flux. 



t(r = i?) =0. 



(2) 



The "hydrodynamic" constraints on the orientation and 
density of long-chain polymer nematics were flrst recog- 
nized by de Gennes [l^, and in the present case, they 
allow us to relate the geometry of ideal chain packings to 
incompressible fluid flow flelds in spherical cavities. 

A divergence-free and azimuthally symmetric configu- 
ration can be described in spherical coordinates as, 



t = mr,9)^ + \7 X 
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where t^{r,9) is any function, and the function '^{r,9) 
must satisfy certain conditions described below. We fo- 
cus first on the class of coaxial spool configurations where 
tcf, ^ and ^ = 0, shown in Fig. [I](a). This configuration 
highlights an underlying element of packing frustration 
in this problem. As t(r) is required to be tangent to the 
sphere at r = i?, the topology of this boundary requires 
the presence point defects, disclinations, of net +2 topo- 
logical charge [201. For axially symmetric configurations, 
these take the form of two -1-1 disclinations at the poles. 
Due to the planar symmetry of the chain orientation in 
coaxial spool, these point defects are spanned by a discli- 
nation line which runs along the axis. Along this defect 
line, the chain flux must vanish in order to maintain a 
single- valued conflguration of t(r), or alternatively, to 
avoid the singular bending energy at the core, leading to 
a physical void at the center of sphere. 

We now consider the second class of configurations, the 
simple solenoid for which t^ = and \I' y ^ 0. ^(r, 9) is the 
well-known Stokes stream function [21[, whose contours 
in the (r, 9) plane describe the orientation of chain flux. 
The components of t(r) are related to the stream function 
via the partial derivatives. 
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Without specifying the details of ^(r, 9) we note that 
boundary conditions place strong constraints on allowed 
conflgurations. From eq. ^ we find that in the limit of 
rsin^ — > 0, we must have that the stream function van- 
ishes as ^ ex r^ sin 9 in order that t(r) remain finite, and 
hence, ^ -^' along the axis. Additionally, the vanishing 
of ir at r = i?, eq. ^, requires that ^ is constant on the 
outer boundary. Since the axis and the boundary meet 
at the poles, these conditions imply that Vl* = along the 
boundary of the sphere and along the axis connecting 
the point disclinations. Configurations of non-vanishing 
^(r, 9) in the volume can be describe in terms of the con- 
tours of constant ^ that cannot intersect and therefore 



form a set of nested (and distorted) tori as shown in Fig. 
[1] (b) . The director lines sweep out the configuration of a 
closed-solenoid, which has the feature that the chain flux 
is non-zero along the polar axis of the spherical volume, 
reducing the disclination line of the spool to two point 
defects in t at the boundary. However, as is clear from 
Fig. [T] (b), this comes at the expense of introducing a 
disclination loop which threads through the interior to 
the solenoid along circular curve where V\l', and hence 
t(r), vanishes. 

The topological requirement of disclination lines in the 
first two families of chain flux begs the question, it is 
possible to achieve non-zero chain flux throughout the 
volume without singularities penetrating into the bulk of 
the polymer nematic? This question is neatly resolved 
by the third class of chain flux configurations, which is a 
superposition of the spool and solenoid with t^ ^ and 
^ 7^ that wc denote as the twisted solenoid. As shown 
in Fig. [1] (c), the lines of t(r) flow along the same class 
of nested-toroidal surfaces of constant \I', but unlike the 
simple solenoid the chain flux also winds around the sym- 
metry axis of the sphere, and unlike the spool t \\ z along 
the axis. Notably, these flux lines have the linked topol- 
ogy of the famed Hopf fibration, a projection of which 
been proposed as model of toroidal condensates of DNA 
and other chiral polymers [22|, [23 , but to our knowledge 
has not previously been studied for its ideal space-filling 
geometry under confinement. For the twisted solenoid, 
t(r) =0 only at point defects at the poles, as spool and 
solenoidal flux configurations "fill in" the voids left empty 
by their counterparts. The chain fiux has a "double- 
twisted" texture around the lines defined by i^ = and 
V* = 0, allowing the defects to "escape to the third 
dimension" along these lines [2^ and crucially, establish- 
ing |t| 7^ throughout the volume. Due to the absence 
of voided disclination lines, the twisted solenoid structure 
generically provides the most uniform flux configurations 
among axially-symmetric, end-free chain packings. 



elastic cost of gradients of the chain flux, both its magni- 
tude and direction, parameterized by the elastic constant, 
K. It is important to note this minimal model includes 
all terms allowed by symmetry to fourth order in t(r) 
and second order in gradients [? ]. 

We analyze the thermodynamics of this non-linear 
model of confined polymer nematics by two methods. 
First, we numerically minimize a finite-difference gener- 
alization of F over all divergence- free, azimuthally sym- 
metric configurations using a standard steepest-descent 
approach on a 20 x 20 grid in the (r, 9) plane. Second, 
guided by the underlying geometrical preferences of the 
model, we construct variational ansatz that allow us to 
study the thermodynamics of the model, albeit approxi- 
mately, in analytical detail. 



A. Coaxial spool 

We first consider the energetics of the spool configu- 
ration in terms of the minimal model of eq. ([S]). To 
a good approximation, this structure is uniform along 
the z direction, and we therefore, study this configura- 
tion in cylindrical coordinates, (p, z, (p), and assume that 
t,f,{v) = tsp(p)- In these coordinates, the free energy den- 
sity has the form 
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While we know of no closed form solution of the Eulcr- 
Lagrange equation for tsp{p), we may construct an ap- 
proximate solution using a well-known approach to non- 
linear field theories, the so-called BPS bound [23|- This 
upper-bound on the true minimum is constructed by 
equating the two positive-definite compressional and elas- 
tic contributions to the free energy, obtaining a first-order 
equation. 



III. THERMODYNAMICS OF CHAIN 
PACKING 

We study the thermodynamics of chain fiux configura- 
tions in terms a minimal free energy that treats the chain 
flux, t(r), as a coarse-grained order parameter field. In 
the spirit of Landau theory, we consider the lowest-order, 
non-trivial expansion of the free energy in terms chain 
fiux and its derivatives subject the divergence free con- 
straint. 



pitl^-tl)±dp{pt^)^0 



(7) 



which has the general solution written in terms of mod- 
ified Bcssel functions t^{p)/tQ = ±A9pln Io{p/X) + 

a{pc)KQ{p/\) . Here, the constant a{pc) = 

—Ii{pc/ X) / Ki{pc/ X) may be adjusted to vary the size of 
a uniformly empty core, with tsp{p < Pc) = 0. However, 
we find that stable spool configurations universally favor 
configurations with pc ^ 0. Consequently, we obtain an 
upper-bound on the spool free energy from. 
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The first term describes thermodynamics of chain pack- 
ing, establishing a constant, preferred value of flux, tg, 
and a corresponding "compressibility" , k, for local devi- 
ations from this state. The second term represents the 



Though the BPS bound approximates the solution the 
Euler-Lagrange equation for isp(p)j it is important to 
note that it rigorously captures the minimal free energy 
configuration in the respective small- and large-p limits. 



This solution for tsp(p) highlights the importance of 
the length scale set by the ratio of compressional and 
elastic moduli, 



/ K \ 1/2 



\Ktl 



(9) 



which increases with the ratio of elastic to compres- 
sive free energy costs. On length scales smaller than 
A, when p ^ A, the elastic cost of filling the discli- 
nation core dominates and the flux at the center van- 
ishes as limp^o *sp/io = p/(2A). On larger length scales, 
when p ^ A, the packing energy dominates and the 
chain flux asymptotically approaches the preferred value, 
limp^oo tsp/to = 1 - A/(2p). 

For the BPS bound on the spool configuration, it is 
easily shown that free energy can be cast purely in terms 
of packing free energy, Fgp = k J dV[tlp{p) — Iq]'^, and 
rescaling all lengths by A we find that the relative stabil- 
ity of the spool configuration is determined by a single 
parameter, R/X. For sufficiently small spherical cavities 
we find that the elastic cost of filling the core outweighs 
the gains in packing free energy and the spool is predicted 
to be thcrmodynamically unfavorable compared to the 
empty sphere with t(r) =0 throughout the volume. For 
the BPS spool ansatz, a first-order transition between 
the empty and coaxial spool configurations is predicted 
at i?ompty/sp = 2.3A, in good agreement with the numer- 
ical calculations which show a transition at 2.0A. 

For larger spheres, R > i?ompty/spj the near-optimal 
flux density at p 3> A stabilizes the coaxial spool. How- 
ever, although the chain density in the spool configura- 
tion approaches uniform filling at large radial distances, 
the slow power-law decay of t^ — t^ ^ p~^ compared to 
linear growth of volume at a given p ultimately leads to 
an appreciable and unavoidable deficiency of the spool 
configuration to fill arbitrarily large spherical volumes 
and thermodynamic cost that grows as ~ i? In R. 



B. Twisted solenoid 

As described in the previous section, the geometry of 
the simple solenoid naturally compliments that of coaxial 
spool in the twisted solenoidal structure, providing the 
most uniform, spacing filling geometry among the three 
chain-packing topologies. For this reason, the coaxial 
spool becomes unstable to the twisted solenoid in the 
limit of highly incompressible chains, or alternatively, in 
the limit of large R/\. To demonstrate this, we construct 
a simple variational ansatz that combines spool configu- 
ration of cq. ([5]) , with a non-zero value of ^ constructed 
to provide a set of nominally, uniformly packed states, 
consistent with boundary conditions. For this we choose 
the following form of the stream function, in spherical 
coordinates 




FIG. 2. Top: four streamline configurations for the twisted- 
solenoid ansatz are shown as A-D. The color scale reflects 
the variation in the magnitude of chain flux, ranging from 
|t(r)| = (blue) in the empty regions to flUed regions 
|t(r)| > to (red). Bottom: the equilibrium value of the 
flux at the sphere center, $, for the minimal model of con- 
fined, end-free polymer nematics. The solid red curve shows 
the results of the variational calculation, while the fiUed- 
points depict results from the finite-difference minimization 
of eq. ([Sjl , both showing a second-order phase transition from 
spool — >■ twisted-solenoidal state. The inset shows plots of 
(|t(r)|)n = (47r)~^/dSl|t(r, n)| corresponding to A-D, the 
profile of flux magnitude at distance from the center, r, aver- 
aged over solid-angle, Q.. 



Here, the cosine dependence on r maintains the condition 
"^{r = R) = and no flux penetrates the boundary, while 



dr[r^i,tw{r)\^r^tl-tl^{r), 



(11) 



-^{r.e) 



r2sin2 6'^twWcos(|^y 



(10) 



is chosen so that, consistent with eq. (|4]), the total 
flux near the core, r ^ R, along the equatorial plane, 
9 = 7r/2, obtains a near uniform value. We analyze the 
following superposition of spool and solenoidal flux con- 



figurations, 



K4/K 

15 10 



ttw(r) = $tsoi(r)+t,p(r), 



(12) 



where tsp(r) is given by eq. ([U and tsoi = {r ^dg'^f — 
dr'^6)/{rsm9). $ is the flux at tlie core of splrerical 
cavity, which serves as an order parameter to describe 
the change of symmetry from the spool state (4> = 0) to 
the twisted solenoid ($ j^O). Fig.[2]shows the streamline 
configuration of this ansatzioi $ G [0, 1], highlighting the 
more uniform flux profile of the twisted-solenoid, which 
reduces the line-like voids of the spool to point-like re- 
gions of depleted chain flux near the poles. 

Inserting the twisted solenoid ansatz into eq. (O we 
have the following free energy for the twisted solenoid, 

i^tw = goiR/X) + 52(i?/A)$' + (74(i?/A)$^ (13) 

a form which is notably consistent with the $ — > 
-$ symmetry. Here, ga{R/\) = Fsp{R/X) 

and g4{R/X) = k/2 J dV\tsoi\'^ positive-definite and 
monotonically-incrcasing functions of R/X, while the 
function 
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encodes both the positive elastic free energy penalty of 
accommodating the solenoidal flux configuration within 
the sphere, as well as the favorable free energy gain 
of back-filling the empty-flux region left behind in the 
spool configuration. The former cost dominates in the 
small- i? limit due to the gradients in t^. and tg (of order 
~ ^to/R) required by the boundary condition, leading 
to cost ~ KR^^. By construction the form of tsoi(r) 
"backfills" the absent fiux underfilled in the spool con- 
figuration, hence, it is straightforward to show that the 
free-energy gain from the more uniform state of chain 
flux also grows as ^ —KtQ^^X^RlnR, always outweigh- 
ing the elastic penalty for large R. Hence, g2{R/X), is 
positive for small R/ and becomes negative above a crit- 
ical sphere size i?sp/tw = 4.7A, indicating a second-order 
phase transition from the spool to the twisted solenoid 
state, 1*1 = ^-32/(2.94) ^ \R- i^sp/tw^/^- Notice that 
in the absence of effects of chain chirality, this phase tran- 
sition represents a spontaneous symmetry breaking and 
is therefore described by the standard mean-field expo- 
nents. 

The dependence of <& on i? predicted by our varia- 
tional theory is shown in Fig. [5J Here, we compare the 
results of our numerical minimization of eq. ([S]), which 
also show a second-order phase transition from a state of 
zero chain flux at the core (spool) to state of increasing 
core flux with increasing container size. We note that the 
numerical results show a phase transition at i? ~ 5.4A, a 
modest shift from the predictions of variational calcula- 
tion despite the rather oversimplified ansatz of both spool 
and solenoidal chain configurations. 




K2tl/K 



FIG. 3. The phase diagram of axially-symmetric, end- 
free configurations of spherical confined polymer nematics in 
terms the reduce sphere radius, R/\, and elastic moduli for 
the non-linear twist and bend costs of eq. ((TSJ . 



C. Twist vs. bend at higher order 

Within the lowest-order, minimal theory of eq. (O, the 
free energy cost of gradients in the chain flux are char- 
acterized by a single elastic constant, K. This is in con- 
trast to the standard Frank elastic energy of oriented 
materials which distinguishes between the relative cost 
of twist deformations (t ■ V x t 7^ 0) and bend deforma- 
tions (t X V X t 7^ 0) [25| . In the polymer- nematic model 
of confined chains, where t(r) is a vector order parame- 
ter, such distinctions will enter the theory as higher-order 
gradient contributions to F. From a microscopic perspec- 
tive of dense chain packing, it reasonable to expect the 
thermodynamic costs of are twist and bend arc quite in- 
dependent. Bend requires intra-chain deformation, while 
twist perturbs inter-chain forces. The underlying physics 
describing these two effects is often distinct, therefore, in 
most cases, not describable by a single phenomenological 
quantity, such as K. We therefore pursue relative contri- 
butions of twist and bend elasticity within in a higher- 
order theory that includes terms allowed to fourth order 
in chain flux and second order in derivatives, leading to 
the following contribution to the minimal model free en- 
ergy, 

5F = ldv{^-^[t-{\/xtf + ^-^[tx{\/xt)]'}, (15) 

where A'2 and K3 arc moduli corresponding to the re- 
spective costs of bend and twist at this order. To study 
the distinct role of twist and bend elasticity in ther- 
modynamics of confined polymer nematics, we analyze 
these higher-order contributions to the free energy based 
on the spool and twisted solenoid ansatz of eqs. ^ 
and ([TU)) . Additionally, we consider a simple solenoidal 
ansatz which is constructed following eqs. ((TU)) - ((T^ . ex- 



cepting for present case of the solcnoidal configuration, 
the azimuthal components of flux vanish (t^ = and 

V-soi = 1/2). 

The higher-order elastic terms have two important ef- 
fects on the gross thermodynamic behavior of the model. 
First, the presence of additional bend costs destabilizes 
the spool configuration with respect to the empty sphere 
configuration, leading to a monotonically increasing de- 
pendence of i?cmpty/sp on K3. For K3 > 6Ktl, the 
suppression of the coaxial spool configuration creates a 
window of stability for the simple solenoid configura- 
tion, which passes from the empty state via a second- 
order phase transition at i? ~ 3.8A. Second, the pres- 
ence of a high twist elastic cost, strongly disfavors the 
twisted solenoid configuration with respect to the un- 
twisted, coaxial spool state. The higher-order elastic cost 
of twist gives rise to a positive contribution to g2{R/X) 
that grows as {K2/2) J dV[t ■ (V x t)]^ - K24R, for 
large R, which is compared to the —kIqR In R dependence 
in the absence of higher-order costs. Hence, the critical 
size at which the spool becomes unstable to the twisted 
solenoid state grows exponentially, as -Rgp/tw ~ e^^''^, in- 
dicating a profound sensitivity of packing phase behavior 
to twist elasticity. 

We show the phase diagram of end-free, confined poly- 
mer nematics in terms of three parameters, R/X, K^^t^/K 
and K2tl/K in Fig [3l The empty, spool and twisted 
solenoid states coexist along a triple line that roughly di- 
vides the thermodynamics into two regimes correspond- 
ing to a high or low ratio of higher-order bend to twist 
costs. A second triple line corresponds the coexistence of 
empty, solenoid and twisted solenoidal states divides the 
phase space into regions of high- and low-bend cost. The 
phase sequence predicted along the line K2 = K3 = 
(discussed above), empty —> spool -^ twisted-solenoid, 
is preserved for sufficiently low K3, although the criti- 
cal values of R/X are shifted as noted above. However, 
when the higher-order costs of bend deformations of the 
chain fiux is sufficiently higher, the coaxial spool phase 
is "squeezed out" by solenoidal phase at intermediate R. 



IV. CONCLUSION 

We conclude with a brief discussion of the parame- 
ters in the minimal model of confined polymer nemat- 
ics, and the relevance to case of dsDNA packing in 
bacteriophage capsids. We crudely estimate the elas- 
tic constant for a dense array of semi-flexible chains 
based on the known cost of bend deflections roughly as 
KtQ « ksTlpd^"^, where ip is the persistence length and 
d is the mean, spacing between chains. Similarly, as- 
suming that the in-plane compressibility of the array 
derives in part from the entropic cost of suppressing 
long- wavelength chain fluctuations (26[ we may estimate 
K^Q ?a ksTlp {d — a)~^/^d~^, a is the chain diameter. 
From these we estimate the fundamental Icngthscale of 
the minimal model as A « ^2/3^1/3^-1/3(2 _. 0i/2)i/3^ 
where </> w (a/^)^ is the volume fraction of the pack- 
ing. For a dsDNA with a = 2 nm and ip = 50 nm in 
densely packed capsid, ^ ~ 0.5 [231, ^^ estimate that 
A ~ 14 nm and a critical capsid size for the transition to 
twist solenoid i?sp/tw — 65 nm. Notwithstanding the sim- 
plicity of the estimate, this length scale is quite compa- 
rable to the typical dimensions of bacterial viruses whose 
radii vary in the range 20 — 100 nm. This correspondence 
as well as the strong 0-dependence of this length suggest 
that the rich thermodynamics of the simplest model of 
single-chain packing plays a necessary role in shaping the 
free energy landscape of genome packing and unpacking. 
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